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PREFACE 


Since the first observation of the ability of a material to absorb 
vibrational energy, engineers and physicists have looked upon this 
property as the possible key to the behavior of materials under dynamic 
Eine conditions. A great amount of work, particularly in the last 
half century, has been done with the purpose in mind of sous 8 
quantitative correlation between these two phenomena. The mere evalua- 
tion of the specific energy absorption property has proven no small task. 
The methods used and units employed in describing this property are many 
and varied so that the uninitiated reader tends to become confused. 

It is the aim of this paper to demonstrate a method of calculating 
the specific energy absorption capacity of a given material. This 
method incorporates the methods of several other experimenters; however, 
it differs in the 5-45 of support of the vibrating member. 

The writer wishes to thank Professors E. K. Gatcombe, R. E. Newton, 
and A. K. Schleicher for their technical assistance and cooperation and 
Machinist Reuel P. Kennicott for his valuable assistance in the fabrica- 
tion and setup of the experimental equipment. 

This work was performed at the U. S. Naval Postgraduate School, 
Monterey, California, during the period January - May 1955 as a require- 


ment for the degree of Master of Science in Mechanical Engineering. 
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CHAPTER 1 


INTRODUCTION 


l. Background. 

The damping capacity of a material is a measure of the ability 
of that material to absorb vibrational energy so that when vibrations 
are induced in a member, the amplitude of vibrations 5522 with 
time without further energy supply, until they become effectively zero. 

For the purpose of this paper, it will suffice to say that the 
energy is dissipated internally within the material by internal 
friction caused by displacements of the crystalline structure under 
load. As yet, the true nature of internal friction is a poínt of 
great conjecture amongst engineers and physicists alike. Many theo- 
ries have been advanced but no firm conclusions as to the microscopic 
behavior of materials in dissipating this vibrational energy have 
been reached. This paper is not intended to attempt an explanation 
of the phenomenon but merely to evaluate its macroscopic effect as it 
concerns the engineer. 

The damping capacity of a material makes itself known in many 
ways --- the visual observation of decrease in amplitude of a vibrating 
spring, the phenomenon of "shaft whirl" where damping causes a lateral 
deflection of a rotating shaft, the increase in temperature of a metal 
member when subjected to cyclic loading of constant amplitude, the 
limited amplitude of vibration of a member when excited by a periodic 
force at its resonant or natural frequency. It is by means of these 
outward manifestations of the damping properties of materials that 


damping measurements are made. 








In general there are six methods that have been used in the past 
to measure damping capacity. These are: 
1. Decay of free torsional vibrations. 


2. Measurement of energy input required to main- 
tain vibrations at constant amplitude. 


3. Measurement of heat generated in a specimen 
vibrated at constant amplitude. 


Lateral deflections of rotating shafts.‏ .با 


5. Forced vibrations at or near the resonant or 
natural frequency of the specimen. 


6. Hysteresis loop in a static stress-strain 
curve. 


Of these methods numbers 1] and 5 have been used to a greater ex- 
tent due to the ease of measurement of the required information. The 
first recorded measurement of damping capacity made in 1906 used the 
method of decay of free torsional vibrations. In recent years 
Foppl (4), Lazan (5), Brophy (1), and Cochardt (2) have used this 
method with good agreement in results. This method is especially 
good if a thin-walled cylinder is used since the stress distribution 
15 then very nearly constant. 

The forced vibrations technique has been used principally by Lazan. 

The units of measurement of damping capacity are quite numerous, 
each experimenter seeming to prefer a different name depending upon 
his method of measurement. This makes reading on the subject rather 
difficult since the reader does not know exactly what the writer is 
referring to when he speaks of the damping capacity of a specific 


material. Some of the more common units are: 





d Viscosity. 

2. Damping capacity. 

Constant of internal friction. 
Hysteretic constant. 

Specific damping capacity. 


Logarithmic decrement. 


Elastic phase constant. 


O =-32 ο νι f$ Ww 


Coefficient. of internal friction. 
In order to remove ambiguity from this paper any references to damp- 
ing capacity will be interpreted as "specific damping capacity" measured 


in inch-pounds per cubic inch of material per cycle of vibration. 


2. Experimental Considerations. 

In any vibrating system there are several ways in which energy may 
be dissipated to produce damping in the system other than the internal 
dissipative ability of the material under investigation. The most evident 
ones are aerodynamic damping due to motion of the system in a fluid medium, 
an energy loss due to energy transfer to the fluid medium by means of 
acoustical radiation, and most important damping will occur in the support 
due to elastic deformation.of the support and/or relative motion of the 
vibrating member with respect to the support. Therefore, any attempt at 
measuring damping capacity must assure that these exterior energy losses 
are of negligible magnitude or that their quantitative effect can be taken 
into account. - 

The method used for the determination of any physical property of a 


material must assure that the material only is being tested and not the 


particular geometric configuration the experimenter may be using. This 





requirement becomes of significance when we consider that the energy con- 
tent or strain energy of a loaded member is a function of the stress level. 
If we are to measure the change in energy we must know how the stress is 
distributed within the loaded member. The method of loading and mode of 
support have a great effect on this stress distribution and therefore 


must be taken into account. 


3. Choice of Experimental Method. 

The method of measurement decided upon by the writer was to deter- 
mine the rate of decay of free transverse vibrations of a cantilever beam. 
In order to minimize the support losses the specimens were made in the 
form of tuning forks. The symmetry of this type of specimen makes it 
possible to consider support losses negligible since deflections at the 
support are extremely small. 

While decreasing the support losses, the tuning fork specimen tends 

to increase the aerodynamic losses due to the relatively broad flat sur- 
faces normal to the direction of motion. The tendency of a tuning fork 

to produce an audible tone when set into vibration also indicates that 
acoustical radiation losses could be high. Since both these losses are 
dependent upon the density of the medium in which the specimen is vibrating, 
the experiment was carried out in an evacuated chamber in order to minimize 
these losses. 

The cantilever beam, a relatively simple mode of support for point 
by point stress analysis, unfortunately has a rather complex stress distri- 
bution. This particular variation is accounted for in the mathematical 
approach to the problem. 

The material tested “in this work was a cold-rolled SAE 1020 steel. 
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4. Summary of Results. 

The only information available for comparison on the damping 
capacity of a 1020 steel was a curve obtained by Lazan. There is no in- 
dication, however, of the condition of this specimen at the time of test. 
Heat treatment and working of the material produce changes in the damping 
capacity as well as the other physical properties so it is impossible to 
make a quantitative comparison between the two curves. The value of the 
damping exponent (slope of the damping capacity versus stress level curve) 
obtained by the writer does, however, fall within the range of values 
given for similar steels. | 

It was Beteruines for the amplitude of vibrations used in this work 
that the losses due to aerodynamic damping and acoustical radiation were 


negligible at atmospheric pressure. 





CHAPTER II 


EXPERIMENTAL SETUP AND PROCEDURE 


1. General. 

The requirement of working in an evacuated chamber to decrease aero- 
dynamic damping and acoustical radiation losses necessitates a more cof- 
plicated setup than would normally be required. 

For the evacuated chamber a 17-inch diameter bell jar was used which 
rested on a steel base plate. The base from a hermetically sealed trans- 
former was set into the base plate to permit bringing in electrical 
leads without destroying the vacuum seal. On the base plate were 
mounted a specimen support, a solenoid for actuating the triggering 


linkage, and the triggering linkage itself. 


2. Specimen Support. 

The specimen support consisted of a steel block two inches square 
and four inches high. It was adjustable in the vertical direction for a 
total travel of one half inch. The top of the block was machined to 
receive the handle of the tuning fork specimen and a top plate, similarly 


machined, was then bolted down to hold the tuning fork securely in place. 


3. Triggering Solenoid. 

The triggering solenoid was a 450 volt A.C. solenoid with a sliding 
core energized from the lighting circuit through a step-up transformer. 
When the solenoid was energized the triggering linkage was actuated to 
rapidly withdraw a spreader from between the tines of the tuning fork 


specimen. 








3. Triggering Linkage. 

The triggering linkage consisted of a single bar pivoted at the lower 
end and attached above the pivot point to the sliding solenoid core. 
The spreader was rigidly attached to the upper end of the bar. (See 


Figures l and 2.) 


h. Specimen. 
The specimen was an accurately machined tuning fork. Machining 
tolerances were such as to give a frequency difference between tines of 


less than .1%. The dimensions of the tines were as follows: 


Length 1,.000" 
Width . 500" 
Thickness 125" 


A lip was machined on the inner surfaces of the tines to permit 


smoother disengagement of the spreader. (See Figure 3.) 


5. Measured Physical Quantity. 

For determining the rate of decay of vibrations, there are two 
physical quantities which can be measured directly. These are the deflec- 
tion of the vibrating member or the surface strains produced by the de- 
flection. The decision was made to measure the surface strains since 
equipment was available for detecting them to’ a high degree of accuracy 
and for transforming them to useable electrical signals. 

The detection of strain was accomplished by mounting a Baldwin SR-4 
type A-1 strain gage on the tuning fork tine. This gage was mounted 
3.500" from the free end of the tine in order to give a maximum signal. 


(Note: Physical dimensions of the gage precluded mounting at the point 
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FIGURE (2) 
Photograph showing specimen support and triggering linkage 
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of maximum strain in the tine.) A gage for temperature compensation was 


mounted on a similar unstressed piece of material. 


6. Instrumentation. 

Accurate instrumentation of any experimental setup is usually one of 
the biggest problems faced by the experimenter. The case at hand was no 
exception. The first attempts were to use a cathode ray oscilloscope to 
present a visual picture of the required information. Since the signals 
from the strain gage are of extremely low voltage, it was necessary to use 
a high gain amplifier between the gage and the oscilloscope. The diffi- 
culty encountered here is that the high amplification required also 
produces considerable noise which distorts the scope trace. Several 
different amplifiers and filtering networks were tried before this method 
was finally abandoned. 

The method finally adopted was to use the Hathaway S-15B oscillograph. 
The Hathaway is a combined electronic, optical, photographic system wherein 
the strain gage signals are electronically amplified and then fed into a 
very sensitive ballistic galvanometer. The movable element of the galva- 
nometer carries a mirror which deflects a beam of light as the movable 
element is deflected in correspondence with the strain gage signals. This 
light beam is focused on a moving sheet of photographic paper which is 
then developed to produce a permanent record. (See Figure 4 for a sample 
record obtained in this manner. ) 


The photograph in Figure 5 shows the complete setup used for this paper. 


7. Procedure. 


The procedure followed for obtaining a record of the decay curve was as 
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FIGURE (5) 
Photograph of complete test setup. 








follows: 


1. Perform all balancing and calibration of detecting 
and recording equipment. 


2. Insert spreader between tines of specimen and 
evacuate bell jar. 


3. With oscillograph photographic paper traveling at 
maximum speed, energize solenoid to remove spreader 
from tines. Maximum speed is used to permit maxi- 
mum discrimination between cycles. 


4. Develop photographic paper. 


8. Analysis of Decay Curve. 

The decay curve was analyzed to determine the logarithmic decrement. 
The logarithmic decrement is defined as the sarita of the ratio of 
two successive amplitudes. Due to the possibility of a very large read- 
ing error occurring in measuring the amplitude of two successive cycles 
of very nearly the same height, an approximating formula was used where 
the decrement was averaged over a number of cycles. The approximating 


& {Am - Am+n 


formula used in this instance was == 7% mmm where n is the 
number of cycles between observations and A n and Amen are the observed 
amplitudes. The error introduced by this approximation is less than 1% 
provided the ratio of amplitudes is less than 1.4. 

The amplitude measurements were made with an optical comparator. The 
reading accuracy of this instrument was .002 inches. 

Logarithmic decrements were calculated over the range of stress from 
approximately 17,000 psi to 3,000 psi to determine the variation of decre- 


ment with stress level. This range of stress was available from a single 


decay trace due to the decrease in stress level as the vibrations decay. 
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CHAPTER III 


CALCULATION OF SPECIFIC DAMPING CAPACITY 


Basically, the method used in this work for calculation of specific 
damping capacity was to determine the specific energy absorption per 
cycle of vibration. 

The information obtained from the decay of free vibrations makes 


possible the calculation of the logarithmic decrement of strain (c) 
Am 


defined as: 8 
e In Am +) 





Since we are interested in an energy change, it is necessary to convert 
this to an energy decrement. 


The logarithmic decrement of strain energy is defined by: 
ሟ ጊን = in Va 


+ | 
Since strain energy is directly proportional to the square of strain, 


Um = In Da y 


the decrement of strain energy becomes: | Ay = In 
Uc vt Ami 


... ain Be = au 
ra! 
With this relation, the energy absorption per cycle becomes: 
-λα 
BU - Un= Un) = σα) 
where Un is the strain energy for the mth cycle of vibration. 


The calculation of Um is accomplished by means of the theorem of 


A ጌ i 
Castigliano: UE Mex) dx 
A, 2ER 


The bending moment M(x) for a freely vibrating beam is a rather complex 
expression since it is a bending moment caused by the inertia forces. 
However, the assumption of a linear variation of bending moment such as 


would be obtained in a cantilever loaded at the free end by a concentra- 
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ted load, will not produce an appreciable error. With this assumption, 


evaluation of Castigliano's theorem becomes: 


2 
A = Eros "ES 
LS h> 
The expression of the change in strain energy now becomes: 


1 

AU = E (CN) 
LSR? 

or, converting to units of stress: 


U= LA ee 
Fh E 





(1) 


Equation (1) now enables us to calculate the total change in strain 
energy of the cantilever. It includes the effects of type of loading 
and support ahd, therefore, does not relate to a specific property of the 
material. 

If we now consider an infinitely ን of material dV, we may 
assume that the energy absorbed in this volume to be 

daus Au dV (2) 
where au is a measure of the specific energy absorption ability of the 
ee (specific damping capacity) per unit volume per cycle. Inte- 
grating, this becomes: 

AU = Jaudv (3) 

B. J. Lazan (5) in nis extensive work on the damping capacity of 
engineering materials has produced consistent results which indicate 


that the specific damping capacity is related to the stress in the 


following manner: 


N | 
Aus KT (4) 


where K and N are constants depending on the partícular material.  Assum- 
ing this equation to be valid, upon substitution in equation (3) we obtain: 
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N 
AU = Ja dv (5) 
In order ‘to evaluate this integral, we employ the bending stress 


equation from elementary strength of materials: 
SES 
ΓΞ με] (6) 
lz 
making the same assumption of linear variation of bending moment as was 


previously made. 


Mac) = May (τ) ! (7) 


ኣ 


This final substitution into equation (5) gives the following 


integral which can be integrated: 


N, N 
NU EE د‎ αν or 


N r4 rh 
Kb M may e N 
e va |» y ho 


Performing the integration and making indicated substitutions (Appendix 1), 
we arrive at the following expression which relates the specific damping 
capacity of the material to the total energy absorption per cycle in the 


vibrating member: 


AU = δαν 
(N 1) (8) 
Solving foryu: | 
AUS (N+) AU (9) AU 


can be calculated from equation (1); however, equation (9) for converting 
to specific damping capacity cannot be solved since it contains the factor 
9698) ),ብ5:5:8 28 ,5(111 én unknown. 

17 








This difficulty can be overcome by referring back to Lazan's equation. 
This equation, when plotted to logarithmic coordinates would produce a 
straight line of slope N. If the derived equation (9) is equivalent to 
Lazan's equation the same result should be obtained from it. 

- Investigation of equation (9) reveals that due to the manner of 
occurrence of the unknown factor N therein, this factor would have no 
effect on the slope of a specific damping capacity versus stress curve 
plotted from πα It is this fact which enables us to determine the 
damping exponent N, by a Simple graphical solution. 

The graphical solution proceeds as follows: An assumed value of 
N - N' is substituted into equation (9) and a fictitious value of 4u = Ou’ 
obtained. su is then plotted against stress to logarithmic coordinates 
and the slope of the resultant straight line measured. The slope of this 
line is then the true value N. Next ce δα intercept is determined 
and corrected by the factor [አላ] to determine the true Au intercept. 
Then the line obtained from the plot of Au' versus stress is transposed 
parallel to itself until it intersects the AU axis at the corrected 
value of AU. The resultant line is now the correct plot of specific 
damping capacity versus stress, and the AU intercept is the constant K in 
Lazan's equation. This method of solution is illustrated in figure (6) 


with N*2. 
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CHAPTER IV 


PRESENTATION OF RESULTS AND CONCLUSIONS 


1. Results. 

The results of this paper, i.e., the damping capacity of a cold- 
rolled SAE 1020 steel, are presented as the curves of Figures (8), (9), 
and (10). Runs #1 and #2 were made in an evacuated chamber in order to 
minimize the losses due aerodynamic damping and acoustical radiation. 
Run #3 was madé at atmospheric pressure for purposes of comparison. The 
equations relating specific damping capacity to stress derived from the 


three curves are: 


#1.4U = 5.25 x 10-5 g 2.30 
#2. AU = 5.28 x 1075 σ-:2:26 
#3. AU = 5.40 x 1075 p 2.29 


These three equations agree within the limits of experimental error. 
The difference between the runs made at low pressure and the runs at 
atmospheric pressure is not large enough to be of significance other than 
the normal scatter and reading error. 

The exponential relationship of Lazan's equation is demonstrated 
quite definitely in all three cases. . 

The point by point evaluation of damping capacity at the various stress 
levels is consistent for the three runs. 

The lack of published datá on the damping capacity of specific materi- 
als make any quantitative comparison of the results of this paper impossible. 
Lazan has published curves of the damping capacity of several steels, one 


of them being & 1020 steel; however, thé condition of the metal at the time 
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Variation of specific damping capacity with stress at a 
pressure of .02 inches of Hg. 
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FIGURE (9) 
Variation of specific damping capacity with stress at a 
pressure of .05 inches of Hg. 
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of test as to heat treatment or forming method is not given. Since we 
would expect changes in damping capacity with heat treatment and forming 
just as in the other physical properties, Lazan's data is not complete 


enough for comparison purposes. 


2. Conclusions. 
The conclusions reached on the basis of these results are: 


l. Lazan's equation relating damping capacity to 
stress level is valid. 


2. The employment of a tuning fork as the vibrating 
.member for determining damping capacities will 
produce consistent results. 

3. For the amplitude of vibrations used in this work, 
the losses due to aerodynamic damping and acoustical 
radiation are negligible. 

3. Sources of Error. 

The most probable source of error in this experimental setup is in 
the analysis of the decay trace to determine the logarithmic decrement. 
The limited speed of the photographic paper in the Hathaway oscillograph 
does not permit adequate sharp definition of the peak of each cycle for 
accurate measurements. The use of a tuning fork designed for a much 


lower natural frequency would remove a large proportion of this source 


of error. 
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APPENDIX I 


EVALUATION OF SPECIFIC DAMPING CAPACITY INTEGRAL 


In Chapter III the expression 


= ει. 
a ff» g^ dedy 


was developed. Integrating this expression and applying the limits we 
obtain: 


u: ib Ma E Ed | 


07 ie. N +1 AUN +!) 


Performing the indicated cancellations this expression becomes: 


AU = Kb ae Baer. 
A (N 


ግ. -N "( IDE 


If this expression is now factored ín the following manner: 





AU = K(bAh)( Mach) s 


it is seen that the term (blh) represents the volume of the rectangular 


N 
η where (Uno BT is v 
liz 


cantilever, and the term K م‎ 
E 9 


& 


represents Lazan's equation for specific damping capacity: 


ኣህ= x 
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APPENDIX II 


Experimental Data and Sample 


Calculations 














39.80 
37.50 
35.45 
33.35 
31.45 
29.60 
27.90 
26.65 
25.10 
23.65 
21.70 
19.65 
17.90 
16.40 
15.15 
14.05 
12.80 
11 5 
10.75 
9.65 

8.45 
7.80 

7.05 


Amen 
37.50 
35.45 
33.35 
31.45 
29.60 
27.90 
26.65 
25.10 
23.65 
21.70 
19.65 
17.90 
16.40 
JE 
14.05 
12.80 
11.75 
10.75 

9.65 

8.45 

7.80 

7.05 

6.40 


RUN #1 


Pressure 


22 
21 
22 
21 
22 
22 
21 
22 
22 
43 
AL, 
43 
43 


LO 
L5 
እሴ 
65 
65 
65 
65 
65 
65 


E 
E 





30 


.02" hg 


. 00270 
.00267 
.00278 
. 00278 
.00275 
. 00269 
. 00218 
. 00272 
, 00276 
. 00190 
. 00225 
› 00217 
. 00204 
. 00195 
. 00188 
. 00207 
. 00195 
. 00137 
.00166 
. 00201 
. 00130 
.00156 


. 00149 





16,550 
15,620 
14,700 
13,800 
13,100 
12,400 
11,750 
11,100 
10,450 
9,550 
8,650 
7,880 
7,220 
6.670 
6,180 
5,630 
5170 
0 وجا 
4,250 
3,720 
3,410 
3,100 


2,820 








41.95 
39.65 
37.10 
34.85 
33.05 
31.50 
30.25 
27.50 
24.45 
22.40 
20.40 
18.65 
a7. 30 
15.80 
14.85 
13.55 
12.55 
12.00 
11:05 
10.30 
9.65 
9.05 

8.65 


39.65 
31510 
34.85 
33.05 
31.50 
30.25 
27.50 
24.45 


22.40 - 


20.40 
18.65 
17.30 
15.80 
14.85 
13.55 
12.55 
12.00 
11.05 
10.30 

9.65 

9.05 

8.65 

8.05 


RUN #2 


Pressure 


22 
21 
21 
21 
20 
20 


Al 


45 


43 


بايا 


43 


hh 


43 


43 


45 





.00256 
.00316 
.00298 
00252 
.00236 
‚00202 
.00232 
.00281, 
.00208 
.00217 
. 00201 
. 00175 
. 00206 


. 00208 


. 00178 


. 00187 
. 00161, 
. 00159 
. 00149 
. 00179 


. 00163 


31 


مهط "03. 


ع ليه 
17,100 


16,150 
15,100 
14,200 
13,450 
12,300 
11,200 
9,960 
9,120 
8,320 
7,600, 
7,050 
6,430 
5,520 
5,120 
1, 500 
1,200 
3,930 
3,680 
3,520 


3,280 


AU 
.0323 


.0343 
«0299 
«0219 
«0185 
«ΟΕ AL 


. 0124 


0074 . 
با 006 . 
0050. 
0038 . 
00318. 


.00311 


.00163 
.00126 
.00107 
› 00068 
. 00097 
. 00075 








39.60 
37.12 
33.20 
21.60 
30.30 
28.20 
27.30 
24.50 
21.00 
19.12 
17.31 
15.80 
13.75 
12.50 
11.18 
10.40 

9.72 

8.72 


37.12 
33.20 
21160 
30. 30 
28.20 
27.30 
24.50 
21.00 
19.12 
17.31 
15.80 
13.75 
12. 50 
11.18 
10.40 
9.72 

8.72 

7.87 


RUN #3 


Pressure = 30.04" hg 


22 


25 


22 


43 


A 


. 00290 
.00301 
400255 
. 00221 
.00253 
.00215 


. 00205 


. 00201 
. 00190 
. 00181 
. 00201, 
.00176 
.00175 
.00160 
«00172 
.00152 
. 00189 


OOLAL 


32 


d ed 
16,000 


14,900 
13,600 
12,500 
11,400 
10,500 
9,700 
8,850 
8,000 
7,250 
6, 500 
5,740 
5,200 
L, 700 
4,350 
1.000 
3,700 


3,380 


0340 
.0298 


.0215 


.0170 


.0168 


.0115 
.0090 
. 0078 
. O061, 
.0052 
.0050 
.0031 
‚0024 
.0021 


.00113 


. 00091 
. 0088 











34. 
.00540 
00521, 
00556 
.00556. 
. 0 
. 8 
. 6 
0051, 
„00510 
.00380 
¿00456 
. 00131 
. 004,08 
.00309 
.00376 

|. 00414 
00390 
.00271, 
. 2 
. 8 
.00312 


. 00296 


- ኋሚ 
e 





. 99161 
«991165 
.9911,5 
. 995 
. 99451 
. 99463 
«99565 
- 99455 
- 994,60 
. 99620 
«99551 
«99565 
. 99590 
«99610 
«99625 
«99585 
«99610 
«92725 
-99670 
99590 
«99690 
› 99700 





.00310 


. 00300 


33 


213 
2 باجا‎ 
216 
190 
172 
154 
138 
123 
109 
91.1 
74.8 
62.2 
521 
h.l 
38.2 
37 
26.7 
22.1, 
18.0 
13.8 
9.6 
7.95 


.0320 
.0281, 
.0260 
.0229 
.0206 
.0179 
.0130 
.0145 
.0128 
.0075 
. 0073 
. 0059 
. 00461, 
.00376 
‚00311 
.00286 
.00226 
„00134 
.00129 
00123 
. 000645 
, 000518 


Sample Tabular Form for Calculation of Specific Damping Capacity (au) 
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